Creation and detection of vortices in polariton condensates 
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We investigate theoretically the creation, persistence and detection of quantized vortices in 
nonequilibrium polariton condensates within a stochastic classical field model. The life time of 
the quantized vortices is shown to increase with the spatial coherence and to depend on the ge- 
ometry of the condensate. The relation with superfluidity in conventional superfluids is discussed. 
Three different ways to measure the vorticity of the polariton condensate are proposed. 



PACS numbers; 03.75.Kk, 05.70.Ln, 71.36.+C. 

The connection between the two intriguing phenom- 
ena of superfluidity and long range spatial coherence 
was conjectured soon after the discovery of superfluid- 
ity. It has proven to be subtle, but is now well estab- 
lished for Bose-Einstein condensates at thermal equilib- 
rium [1]. The recent experimental observation of the 
spontaneous build up of long range spatial coherence in 
exciton-polariton gases [2] as well as a bimodal momen- 
tum distribution ^,t2i,t4i] and long temporal coherence 5], 
has induced the interest in the phenomenon of Bose Ein- 
stein condensation out of thermal equilibrium. Due to 
the short (ps) life time of the polaritons, they have to 
be continuously injected in the microcavity by means of 
relaxation from an excitonic reservoir that is excited by 
a laser. The resulting steady state is not a thermal equi- 
librium one, yet it still features the main characteristic 
of an equilibrium Bose-Einstein phase transition: short 
spatial coherence below a critical density, long range co- 
herence above. A natural question is whether the phase 
with long range coherence manifests superfluid proper- 
ties. The identiflcation of a superfluid regime out of 
equilibrium has been addressed for a resonantly driven 
polariton gas, where the coherence is imprinted by the ex- 
citation laser [6, 7]. Under this excitation condition, su- 
perfluidity manifests itself as a reduction in the Rayleigh 
scattering when the polariton gas is injected at finite mo- 
mentum. In this Letter, we will address the issue of su- 
perfluidity of a nonresonantly excited polariton gas where 
the long range order is spontaneously formed and super- 
flows are not directly forced by the excitation laser. 

Phenomenologically, superfluidity can be loosely de- 
scribed as 'flow without friction', referring to the re- 
markable property that the superfluid does not neces- 
sarily move along with the container that holds it. This 
idea of 'flow without friction' does not apply to driven- 
dissipative systems: When energy is fed into the system, 
it does not come as a surprise that the fluid is in mo- 
tion with respect to its container. This is why it was 
not possible to interpret the observation of the quantized 
vortices in polariton condensates as a proof of super- 
fluidity. The appearance of vortices in the steady state of 
a polariton condensate was rather interpreted as a result 
of the flow driven by the relaxation of high energy ex- 



citons into the lower polariton branch that is subject to 
an irregular potential, a phenomenon that could equally 
well be exhibited by a classical system. The same holds 
for the spontaneously formed vortex lattice that was pre- 
dicted in the theoretical work of Ref . [0] . 

A slightly more general view on a superfluid is that it 
supports many metastable flow patterns under the same 
boundary conditions. Such behavior is not displayed by 
classical systems: the velocity of water in a pipe with 
a given pressure difference between both ends converges 
quickly to its unique steady state. This definition of su- 
perfluidity as multiple metastable flow patterns is equally 
meaningful in nonequilibrium systems. We will propose 
an experiment to investigate the superfluidity of polari- 
ton condensates in this sense and make its theoretical 
analysis within the frame work of the classical field model 
introduced in Ref. [lo| . In particular, we will investigate 
the evolution of an additional vortex introduced in a vor- 
tex free polariton condensate. 

A mechanical creation of vortices, in analogy with su- 
perfluid ''He and ultracold atomic gases [ul, is proba- 
bly not practical, because it would require mechanical ro- 
tation frequencies between 1 GHz and 1 THz, but thanks 
to the half light nature of polaritons, a vortex can be cre- 
ated with a pulsed 17| resonant laser beam that contains 
one unit of angular momentum per photon. Because it 
is at present technologically impossible to perform single 
shot measurements, direct visualization of a moving vor- 
tex in the density or phase profile of the condensate can- 
not be achieved. We will propose instead three different 
operational methods to detect the condensate vorticity 
in multiple shot measurements. 

Our classical field model describes the coupled dynam- 
ics of a classical incoherent exciton reservoir to the quan- 
tum dynamics in the lower polariton branch. The cou- 
pled equations of motion of the reservoirs describe the 
nonresonant excitation with intensity P, and the scat- 
tering into the lower polariton region 
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The coupHng terms between the reservoir and lower po- 
lariton branch read explicitely 

q,x' 

X i?i„,out(eq)V'(x') (2) 

The reservoir density is normalized so that threshold 
is reached when it is unity. The lower polariton dynamics 
is described by the stochastic classical field equation 
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The noise term dW is a complex Gaussian stochastic vari- 
able with the correlation functions: 
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where nf„^^^ 



[7?.in,out + (7?.in,out)'^]/2 are the sym- 
metrized in and out-scattering kernels. According to 
Eq. ([3]), the number of excitons that scatter from the 
reservoir into the lower polariton classical field equals 
^(-0'''V')|res = 2Re[V'*(7?.in-7?.out)'0]- The parameter a in 
Eq. H]) quantifies the strength of the coupled to the lower 
polariton field V'! the last term in Eq. ([1]) comes from 
the fact that the classical field in Eq. ^ contains 1/2 
zero point fluctuations per mode, because it samples the 
the Wigner quasi- probability distribution Pw, that is re- 
lated to the observable density as / Pw{i')'^P*4' = + 
'0V'^)/2. In the simulations, we have used the following 
parameters: 7 = 0.5 meV, rriLp/h = 1 /im^^meV^^, 
g/h = 0.005 meV^Ltm^, Rout{E) = (3meV), Rin{E) = 
eyipi- E/kBTR)[-f + 3meV], with keTR = 2 meV. Nu- 
merical results are obtained by averaging over 100 real- 
izations of the noise. 

Vortices are typically more stable in a toroidal geome- 
try as compared to a simply connected one. We therefore 
start the analysis of the vortex life time for a polariton 
gas that is excited by a ring shaped excitation laser. The 
right hand panels in Fig. [1] show the simulated evolu- 
tion of the angular momentum and excess density when 
a Gauss-Laguerre beam with one unit of angular momen- 
tum is applied to the polariton gas for three different con- 
densate densities (increasing from top to bottom). The 
corresponding spatial coherence is shown in the left hand 
panels. The excess density (dashed line) evolves almost 
identically in the three cases [3], but the time depen- 
dence of the angular momentum is dramatically differ- 
ent. At low densities, the angular momentum quickly 
decays together with the excess density introduced by 



the probe pulse, where at the highest density, the angu- 
lar momentum lasts up to the end of our simulation time 
window (120 ps). The insets in the left column show the 
spatial coherence at the final time of our simulation win- 
dow. In the upper and central panels, where the vortex 
disappears at the final time of the simulations, an im- 
portant decrease in the absolute value of the coherence 
is observed. 

The simultaneous increase in the coherence and the 
vortex life time is in analogy with conventional superflu- 
ids. The argument for the difference in life time of the 
angular momentum is that when the phase of the conden- 
sate is coherent over large distances, a local distortion of 
the phase is not possible and vorticity can only disappear 
if the vortex core moves out of the polariton gas. The 
time it takes for a vortex to move out depends on the ge- 
ometry of the condensate. In polariton condensates, an 
important phenomenon that determines the life time of 
the superflow are the currents that arise due to the inter- 
play between the nonresonant excitation and the external 
potential [13| . Vortices are dragged by these currents and 
in order to achieve a long life time of the vortex it is ben- 
eficial if these currents are pointed toward the center of 
the vortex. Under excitation with a ring shaped excita- 
tion laser of Fig.[Tl the blue shift causes the polaritons to 
flow toward the center, which hinders the vortex to leave 
the system, because it has to move upstream. Fig. [2] 
shows the same simulations with a regular top hat ex- 
citation laser. The polariton-polariton interactions now 
cause an outward instead of inward flow and the angular 
momentum in the phase coherent regime decays much 
quicker. 

The phase rigidity of the polariton condensate becomes 
also clear in the response of the angular momentum for 
varying intensity of the Gauss-Laguerre pulse. At low 
intensities, the pulse is not able to create a vortex in the 
condensate so that the transfer of angular momentum is 
poor and decays quickly. Only when a certain threshold 
intensity is passed (corresponding in our simulations to 
an excess density of about half the condensate density) a 
vortex is formed and the angular momentum persists for 
long times. This behavior is illustrated in the inset in the 
lower right panel of Fig. [U where the response to a low 
intensity pulse is shown. In the normal, incoherent phase, 
we have observed no qualitative change in the transfer of 
angular momentum: it always decays on the same time 
scale as the excess density introduced by the pulse. 

We now turn to the question how the metastable super- 
flow can be detected experimentally. Ideally, an experi- 
mental measurement should be time resolved, so to ob- 
serve the predicted increase in life time of the metastable 
superflow when the condensate density is increased. 

A first technique is the one that was used previously 
to detect vortices induced by the disordered potential 
acting on the polaritons: A vortex leads to a fork dislo- 
cation in the interference pattern of the condensate with 
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FIG. 1: The left column show the first order spatial 
coherence between inversion symmetric points g^^\x) = 
(i/)^ (— x)i/)(x))/ ^n(x)n(— x) for values of the density increas- 
ing from left to right (n = 1.3, 5.7 and 32 /im~^ respectively), 
before the arrival of the Gauss-Laguerre pulse; the insets show 
the values at the latest simulated times. The right column 
shows the temporal evolution of the total polariton density 
(dashed line, normalized to its maximal value) and angular 
momentum (thick) after a Gauss-Laguerre probe pulse is has 
perturbed the polariton condensates. The thin line shows the 
observable A that we propose to use for a time resolved mea- 
surement of the vorticity. The inset in the lower right figure 
shows the response of the condensate to a probe pulse that 
is too weak to create a vortex. The polariton condensate is 
created with a ring shaped excitation laser (inner radius 7/im; 
outer radius 20fim). 



its inverted copy [8|]. In order to observe the increased 
life time of the superflow when the condensate density 
is increased, this interferogram should be measured with 
time resolution. This is possible with the use of a streak 
camera. Still, an average over many runs of the experi- 
ment should be recorded in order to have a good signal to 
noise ratio. This poses some difficulties when the vortex 
undergoes a random motion that is different from shot 
to shot. In that case, no dislocation will be visible in an 
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FIG. 2: The same as Fig. 1, but for a condensate excited 
with a top-hat excitation laser with a 20/im radius (maximum 
density n = 48 /im~^). 



interferogram that is averaged over many realizations. 

A first possibility to do a time resolved measurement 
is based on the change of sign of g^^^ (x, — x) (compare 
the left hand plots in Fig. [1] with their insets). In a 
time-resolved interferogram, such as the one recorded in 
Ref. [3] , this would be visible as an interchange between 
light and dark fringes as a function of time. The random 
motion of the vortex leads to a reduced contrast of the 
fringes, but fringes are not washed out (the coherence 
between two points at opposite positions of the vortex in 
the right panel of Fig. [T] is about 0.4). 

A second possibility is a direct measurement of the 
time dependence of the angular momentum. We discuss 
here how this can be achieved with interferometry. The 
operator that rotates a state over an angle 6 can be writ- 
ten as Rg = eKp(i9Lz). Letting the polariton condensate 
interfere with a copy that is rotated by an angle 9 and 
phase shifted by a leads to a integrated intensity, nor- 
malized to the total number of particles N 
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For small rotation angles and a w 7r/4, the argument of 
the cosine can be expanded. It is then easy to see that 
the contour line of the interferogram I{a, 6) = 1 goes 
to {Lz)9 = tt/2 — a. If the fluctuations of the angular 
momentum are small, the contour lines of the interfero- 
gram are straight lines and the angular momentum can 
be determined from a large region in the (a, 9) plane. 
Fig. [3] shows simulated interferograms. Left and right 
panels are before and after the Gauss-Laguerre pulse re- 
spectively. In the right hand panel, the vorticity of the 
condensate results in a tilting of the interferogram in the 
(a, 9) plane. Note the decrease in the contrast of the 
interferogram at finite rotation angles 9. This is a conse- 
quence of the imperfect spatial coherence. 

An alternative scheme that could be used to detect 
the presence of a vortex is based on the fact that a vor- 
tex induces a flow of polaritons v cx le/r. Therefore, 
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FIG. 3: Interferogram as a function of the phase difference a 
and rotation angle Q betweeen the copies of the condensate 
that are interfered, before (t=0 ps, left hand panel) and after 
(t=100 ps, right hand panel) the arrival of the Gauss-Laguerre 
pulse. The black lines show the contour lines J(a, = 1; the 
white line shows a = 7r/2 + {Lz)0, where (L^) is taken from 
Fig. [T] The presence of angular momentum is signalled by 
the tilting of the interference fringes in the (a, 9) plane. 

measuring the velocity (momentum) of polaritons on a 
path around the vortex can reveal its presence. The pro- 
posal to measure this velocity is sketched in Fig. [D It 
consists of measuring the momentum distribution of a 
slice of polaritons selected in real space. The presence 
of vortex will be signaled by a shift in the momentum 
distribution along the long directions of the slits. The 
observable that corresponds best to the angular momen- 
tum is A = {ks — few — kw + ks). It is plotted in Figs.[T] 
and [5] together with the angular momentum L, showing 
a good correspondence between A and L. Note that this 
technique works better in the case of a regular top hat 
excitation laser, where the vortex sits in a high density 
region. In the ring excitation geometry, this detection 
technique is more difficult, because no region with simul- 
taneously a sizable density and velocity can be found. 
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FIG. 4: Sketch of the scheme to detect the extra velocity in 
the polariton condensate due to the presence of a vortex: The 
momentum distribution from a slice selected in real space is 
shifted in the presence of a vortex (full line) with respect to 
the momentum distribution in its absence (dashed line). 

In conclusion, we have theoretically investigated the 
possibility to create and detect quantized vortices in po- 



lariton condensates by resonantly exciting them with a 
Gauss-Laguerre laser beam. Only when the system ex- 
hibits a large degree of long range order, the life time 
of the angular momentum substantially exceeds the life 
time of the excess density introduced by the Gauss- 
Laguerre beam. We have proposed and analyzed sev- 
eral schemes to experimentally detect the presence of a 
quantized vortex. 
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